A Model for Superconductivity in Ferromagnetic ZrZn 2 
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This article proposes that superconductivity in the ferromagnetic state of ZrZn2 is stabilized by an 
exchange-type interaction between the magnetic moments of triplet-state Cooper pairs and the fer- 
romagnetic magnetization density. This explains why superconductivity occurs in the ferromagnetic 
state only, and why it persists deep into the ferromagnetic state. The model of this article also yields 
a particular order parameter symmetry, which is a prediction that can be checked experimentally. 
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Recently, superconducting states have been found to 
coexist with ferromagnetism in the materials UGe2 [fij 0, 
§, ZrZn 2 {§, and URhGe §. The initial discovery in 
UGe2 was motivated by the idea that parallel-spin (and 
hence spin-triplet-state) Cooper pairs would be favored 
in a metallic state close to the border of ferromagnetism. 
The proximity of a ferromagnetic state would give rise to 
relatively strong ferromagnetic fluctuations which would 
promote spin-triplet pairing. 

A sketch of the phase diagram as measured in Ref. || 
for ZrZn2 is shown in Fig. R|. As noted by the authors 
of Ref. [|, one of the most intriguing and perhaps sur- 
prising features of superconductivity in ZrZn 2 (as well 
as that occurring in UGe2) is that it occurs only in the 
ferromagnetic state. The reason for the surprise is that 
previous theoretical work had not anticipated that su- 
perconductivity could occur in the ferromagnetic phase, 
unless at the very least it was also stable in the para- 
magnetic phase. The possibility that superconductivity 
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FIG. 1: Phase diagram showing the ferromagnetic (T/) 
and superconducting (T s ) transition temperatures in ZrZn2 
as functions of pressure, as derived from the model of this 
article and as determined by experiment. Q For clarity, the 
temperature scale for the superconducting phase transition 
has been multiplied by a factor of approximately 10 relative 
to that for the ferromagnetic phase transition as in Ref. ^ 
Note that the qualitative behaviors of both Tf and T s for 
UGe2 (sketch of data from Ref. ^) are quite different from 
those for ZrZn2. 



might appear only in the ferromagnetic phase does not 
seem to have been considered before the recent experi- 
mental discoveries. For example, a very early article j^] 
had noted that the presence of the large internal mag- 
netic induction in a ferromagnet would suppress super- 
conductivity. Also, another early theoretical article [Q, 
which demonstrated how spin fluctuations can give rise 
to p-wave superconductivity, found that the supercon- 
ductivity occurs in both the ferromagnetic and paramag- 
netic phases close to the ferromagnetic quantum critical 
point (see Fig. ||). Other examples which find supercon- 
ductivity on the paramagnetic side of a ferromagnetic 
quantum critical point include Refs. ||, ||, [To[ [ll]. Very 
recently it has been argued (l2) that the critical tem- 
perature for spin-triplet p-wave superconductivity medi- 
ated by spin fluctuations is generically much higher in 
the Heisenberg ferromagnetic phase than in the param- 
agnetic phase, due to the coupling of the magnons to the 
longitudinal spin susceptibility, and this result is qual- 
itatively in agreement with the superconducting phase 
diagram for UGe2 (see Fig. |J) . Another line of argument 
]l3| is that the pairing symmetry realized in UGe2 must 
be a nonunitary spin-triplet pairing similar to that real- 
ized [pj in the Ai-phase of superfluid 3 He because such 
states are free from the Pauli limit and can survive in a 
huge internal magnetic field. In addition, the supercon- 
ducting order-parameter symmetry in the ferromagnetic 
phase of UGe2 has been studied in terms of the mag- 
netic point group symmetry of the ferromagnetic phase. 
The ideas introduced below have some overlap with these 
latter |p| 14, |l5| ideas. Finally, we note an article that 
has shown theoretically that coexisting superconductiv- 
ity and ferromagnetism can occur for the case where the 
same band electrons produce both phenomena E3|. 

This article describes a phenomenological model that 
gives a good description of superconductivity in ferro- 
magnetic ZrZn2 (although not in ferromagnetic UGe2). 
In particular, the model gives a natural explanation of 
the fact that superconductivity occurs in the ferromag- 
netic but not in the paramagnetic phase. The basic idea 
is that in the superconducting state the Cooper pairs can 
have magnetic moments — see Ref. |l7]. In the presence 
of a ferromagnetic magnetization density in the sample 
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FIG. 2: Schematic reproduction of the results of Ref. ^ (solid 
line) and Ref. [l^ (dashed line) showing theoretical calcula- 
tions of the p-wave superconducting transition temperature 
versus pressure P. (Here P represents any parameter charac- 
terizing the distance from the quantum critical point.) The 
result of Ref. |^ and others (see text) were responsible for the 
idea that superconductivity in a ferromagnetic state would 
also be accompanied by superconductivity in the neighboring 
paramagnetic state. The very recent result of Ref. fL2j is qual- 
itatively similar to the phase diagram determined for UGe2 
(which has a small pocket of superconductivity close to P c ) 
but not to that for ZrZn2 (where superconductivity occurs at 
all P between P = and P = P c ). 



the magnetic moments of the Cooper pairs can inter- 
act with this ferromagnetic magnetization density via an 
interaction having the form of an exchange interaction. 
The Cooper pair magnetization density chooses a direc- 
tion that makes this "exchange" energy negative, and 
this is the mechanism that makes the superconducting 
state more stable in the ferromagnetic state than in the 
paramagnetic state. As will shown below, in order to give 
rise to superconductivity in the ferromagnetic state but 
not in the paramagnetic state, the exchange coupling just 
described must be greater than a certain critical value. 

The ferromagnetic state will be modelled using the 
Landau free energy Q 



F, 



1 



< f = a' f [T-T f (P)]M* + -0 f M 



(1) 



Here the ferromagnetic transition temperature Tf is as- 
sumed to depend on the pressure P. Expanding Tf(P) 
in a Taylor series about the point P c at which it goes to 
zero, and keeping only the first nonvanishing term, yields 
T f (P) = T' f {P c ~ P). This linear dependence of T f on 
P agrees well with the experimentally measured pressure 
dependence for ZrZn2 [0, shown in Fig. (D. From Eq. (|l|) 
one finds M = {a' f / f3 f y/ 2 (T f {P) - T) 1 / . 

For cubic ferromagnets (such as ZrZn2) the only two 
possibilities for the easy direction of the ferromagnetic 
magnetization density are the a [100] or a [111] direction 
p8|. Although in the absence of ferromagnetism, the 



C15 Laves phase structure of ZrZn2 has cubic (m3m) 
point group symmetry, the point group symmetry in the 
presence of ferromagnetism is reduced to the magnetic 
point group D^i^C^) (4/mm'm') symmetry for the fer- 
romagnetic magnetization density in the [100] direction, 
or T>3d{C3i) (3m') for the ferromagnetic magnetization 
density in the [111] direction. Since all of the irreducible 
representations of C^h and are one dimensional, it 
is expected that the transition to superconductivity in 
the presence of ferromagnetism can be described by a 
one-component order parameter. 

It is of interest to investigate how this one-component 
order parameter describing superconductivity in the pres- 
ence of ferromagnetism might be related to order pa- 
rameters appropriate to the description of supercon- 
ductivity in cubic ZrZn2 in the absence of ferromag- 
netism. An advantage of treating the paramagnetic, 
non-superconducting state as the reference state is that 
an explicit dependence of the parameters describing the 
superconductivity on M is obtained (see below). Be- 
cause of the large value of the exchange field compared 
to the superconducting critical temperature, all spin- 
singlet states of Cooper pairs are excluded. Thus, as- 
suming spin-triplet pairing, consider the representation 
Flu of the group for which the order parameter is 
the three-component quantity i/> = (tp x ,ip y ,tp z ) whose 
components transform under rotations like those of a 
three-dimensional polar vector |l9| (the F^ u representa- 
tion gives the same model). We use a strong spin-orbit 
coupling scheme, in which rotations transform both spin 
and orbital degrees of freedom. Also, the time-reversed 
state corresponding to xjj is tb R — Now de- 

fine the vector product 



S = iip* x 



(2) 



Because this quantity transforms like a magnetization 
density under the operations of and time reversal, 
it will be interpreted (to within a constant factor) as a 
magnetization density associated with the Cooper pairs. 
It should be noted that, at the phenomenological level of 
this article, in the strong spin-orbit coupling scheme, the 
spin and orbital magnetization density of Cooper pairs 
can not be distinguished. 

Now consider the following terms of an expansion of 
the Ginzburg-Landau (GL) free energy in powers of the 
components of the order parameters M and if) 



F 



atp* ■ - 4-kJM ■ S. 



(3) 



Only the terms quadratic in the superconducting order 
parameter and consistent with the cubic symmetry and 
time-reversal invariance have been included here, since 
these are all that are necessary (together with the gra- 
dient terms) to find the upper critical field for super- 
conductivity. Furthermore, terms up to linear order in 
M have been included. Note that the last term in this 
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equation has the form of an exchange interaction between 
the ferromagnetic magnetization density and the Cooper- 
pair magnetization density. If the exchange parameter 
J is positive, the formation of a superconducting state 
in which the Cooper-pair magnetization density is paral- 
lel to the ferromagnetic magnetization density is favored. 
Also, if the ferromagnetic magnetization density M is ro- 
tated (by an applied magnetic field) this exchange mech- 
anism for stabilizing the superconductivity is still appli- 
cable, and the orientation of the Cooper-pair magneti- 
zation density S will follow that of the ferromagnetism. 
(If J < 0, an equivalent model is obtained in which a 
Cooper-pair magnetization density antiparallel to the fer- 
romagnetic magnetization density is favored.) The free 
energy of Eq. (||) is reminiscent of that employed in the 
description of the Ai-phase of 3 Hc fl4| . 

Now call the direction of the nonzero ferromagnetic 
magnetization density the z direction (which, as noted 
above, can be either a [100] or a [111] direction for a 
cubic ferromagnetic). In addition to the exchange field 
coupled with the spin of electrons, the magnetization cre- 
ates an internal magnetic induction which interacts with 
the electron charge. Thus, the superconductor should be 
in the mixed state even in the absence of an external 
magnetic field, and, in order to calculate the transition 
temperature, one has to take into account the gradient 
terms in the GL free energy, in addition to the uniform 
terms given by Eq. (||): 

F s = F S:0 + ^(A^J'fA^) + if 2 [(A^)*(D^ 3 ) 
+ {D l ^)*(D 3 yj l )} + X 3 (AV^)*(A^)- (4) 

Here a = a'(T — To), and To is the superconducting 
transition temperature in the absence of the exchange 
interaction of Cooper pairs with the ferromagnetic mag- 
netization (i.e. at J = 0), which is assumed to be posi- 
tive. The gradient part contains terms which are invari- 
ant under rotations from the cubic group |2(| , with Di = 
— ih(d/dxi) + (2\e\/c)Ai, and curL4 = B. The magnetic 
induction is given by B = H ext + AitM ', where H ext 
is the external magnetic field directed along z, and the 
magnetization density is M = Mo + (// — \)H ext / (Air). 
The long cylinder geometry, with the z axis along the 
axis of the cylinder, has been assumed. 

Using a variational approach pl| ] to minimize the free 
energy (||), we calculate the superconducting critical tem- 
perature as a function of the magnetization density and 
external field, which takes the simple form 



T C (M) = T 



MJ - Jc) 



A I 



(5) 



at H ext — 0. The quantity J c describes the suppression 
of the critical temperature due to orbital effects. It takes 
different values for M || [100]: J c = {\e\/hc){2K 1 + 2K 2 + 
K 3 ), and for M || [111]: J c = (\e\/hc)(2K 1 + 2K 2 + 
2K3/3). Another result of our calculation is that the only 



component of the order parameter which is non-zero at 
T = T C (M) - is 0- = (Vfc ~ i^ y )/V2 (or V+ = (^ ± 
iip y )/V2 for J < 0). It is this quantity that describes the 
formation of a superconducting state with its Cooper- 
pair magnetization density parallel to the ferromagnetic 
magnetization density. Finally, note from Eq. (^|) that, in 
order for T c to be enhanced in the ferromagnetic phase 
relative to its value To in the paramagnetic phase, the 
exchange parameter J must be greater than J c . In the 
weak-coupling theory, J is proportional to A'(ef), the 
derivative of the single-particle density of states (DoS) 
at the Fermi level p7| . The smallness of this quantity 
in 3 He explains the narrow region of existence of the Ai- 
phase. In the case of ZrZn 2 , however, where the DoS 
is extremely sharply peaked near the Fermi energy [p2[ , 
N'(cf) could be very large, but estimating J in terms of 
N'(6f) is probably too simplistic. 

In order to confirm tp_ as a possible order parameter 
describing the formation of superconductivity in the fer- 
romagnetic state, it should be checked that it transforms 
as a basis vector of some irreducible representation of 
the magnetic symmetry group of the ferromagnet [ jl5| . 
Suppose that the ferromagnetic magnetization density is 
along the [100] direction. Then the magnetic symmetry 
group is B 4h (C4 h ) = C 4h + (RC 2x )C 4h . (Here R is the 
time-reversal transformation.) In this case, ip- trans- 
forms like one of the complex irreducible representations 
( 1 E or 2 E) of C4/1. Furthermore, although there is no 
time reversal operation in this magnetic group, the op- 
erator RC 2x has the effect of replacing ip_ by its com- 
plex conjugate. Hence ip_ is a possible order param- 
eter. A similar analysis can be performed if the fer- 
romagnetic magnetization density is along the [111] di- 
rection, when the magnetic point group is D^iCsi) = 
C31 + (RC 2x )C3i. Here too the order parameter trans- 
forms like one of the complex representations ( 1 E or 2 E) 
of the relevant point group (C3i). It should be noted that 
this predicted symmetry of the superconducting state can 
be verified by experimental measurement (e.g. sec [p3[ ). 

The pressure dependence of the critical temperature 
T s of the transition to the superconducting state can be 
found from Eq. (JsJ) to be given by the solution of 

T s = T + T* 1/2 [Tf(P)-T s } 1/2 



where 



a' J 



(6) 



(7) 



By assumption, the exchange enhancement results in a 
superconducting transition temperature T s much greater 
than the superconducting transition temperature To in 
the paramagnetic state. Furthermore, except for P very 
close to P c , T s <§; Tf(P). Under these conditions the 
pressure dependence of T s is given by the formula 



T S {P) = T S {0)(1-P/P C ) 



1/2 



(8) 
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When P gets very close to P c and Tf(P) becomes very 
small this equation is no longer valid. In this extreme 
circumstance, if one takes T = and Tf (P) -C T* , one 
finds 



T s (P)=T f (P)[l-T f (P)/T* + ...} 



(9) 



which shows that, for P very close to P c and To = 
0, T S (P) approaches Tf{P), and is never greater than 
Tf (P) . Eq. (||) , together with the equation for the ferro- 
magnetic transition temperature 



Tf(P)=T f (0)(l-P/P c ) 



(10) 



arrived at above, have been used to plot the phase di- 
agram of Fig. (Q), which shows a remarkable similarity 
to the experimentally determined 0] phase diagram for 
ZrZn2. 

Finally, at a given pressure, we find a temperature de- 
pendence of the upper critical field of 



H c2 (T, P) = H c2 (0, P) [1 - T/T s (P)] 



(11) 



This is in reasonable agreement with the experimental 
result (see Fig. 3 of Ref. |4j , which however has somewhat 
more curvature than the linear temperature dependence 
shown here. The lack of curvature in the result of Eq. 
( pi] ) results from the linear dependence of M z on H ext in 
our relation M z = M z q + (p— l)H ext / (4tt). 

We conclude that the proposed mechanism of stabiliz- 
ing superconductivity in a ferromagnet (by an exchange 
type of interaction between the magnetization density of 
the Cooper pairs and the ferromagnetic magnetization 
density) gives an excellent qualitative description of the 
phase diagram determined experimentally for ZrZn2. In 
particular, it explains in a natural way the fact that the 
superconductivity occurs in the ferromagnetic phase, but 
not in the paramagnetic phase. For this mechanism to 
work, the exchange interaction parameter must have a 
magnitude larger than a certain critical value. A further 
experimental test of our model would be the determina- 
tion of the order parameter symmetry. (A prediction of 
our model is that the order parameter transforms like 
one of the complex representations of the relevant point 
group.) 

It should be mentioned that the spin-fluctuation mech- 
anism studied in Ref. |l2[ can provide an alternative ex- 
planation of growing T s in the ferromagnetic state, given 
that the magnetization in ZrZn2 does not reach satura- 
tion. To what extent the fluctuation effects discussed in 
Ref. |E| are essential compared to the mean field inter- 
actions studied in this article, is in our view still an open 
question, and their relative contributions can be differ- 
ent in different materials. For example, it seems that the 
phase diagram of UGe 2 (see Fig. 0) can be satisfactorily 
explained by the spin-fluctuation theory, and the appar- 
ent absence of contributions from the exchange interac- 



tion of our work could be explained by the magnitude of 
the exchange parameter being less than its critical value. 
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